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We calculate the screening corrections to the electron capture rates in dense 
stars by the relativistically degenerate electron liquid. In order to calculate the 
screening corrections we adopt the linear response theory which is widely used 
in the field of solid state physics and liquid metal physics. In particular, we 
use the longitudinal dielectric function for the relativistically degenerate electron 
liquid derived by Jancovici. We calculate the screening potential at the position 
of the nucleus. By using this screening potential one can calculate the screening 
corrections to the electron capture rates. We will present accurate analytic fitting 
formulae which summarize our numerical results. These fitting formulae will 
facilitate the application of the present results. The screening corrections to the 
electron capture rates are typically a few percent. 
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1. INTRODUCTION 

Since the pioneering works of Fuller, Fowler, & Newman (1980, 1982a, 1982b, 1985), 
calculations of stellar weak-interaction rates have entered an era of precision science. More 
recently an authoritative work of Langanke & Martinez-Pinedo (2000) on this subject ap- 
peared. 

Since the presupernova stellar evolution and stellar nucleosynthesis critically depend on 
the details of the stellar weak-interaction rates (e.g., Wanajo et al. 2002), it is extremely 
important to calculate accurately the screening corrections to the electron capture rates in 
dense stars by the relativistically degenerate electron liquid. 

This problem has been already addressed by many authors (Couch & Loumos 1974; 
Takahashi, El Eid, & Hillebrandt 1978; Gutierrez et al. 1996; Luo k Peng 1996; Bravo & 
Garcia-Senz 1999). The plasma effects on the chemical potential of the nucleus and hence on 
the threshold energy for the electron capture, in particular, have been discussed by Couch 
& Loumos (1974), Gutierrez et al. (1996), as well as by Bravo & Garcia-Senz (1999). 

In this paper we will address ourselves to the calculation of the effective potential energy 
felt by the relativistically degenerate electron. We will use the linear response theory in order 
to calculate the screening potential caused by the relativistically degenerate electron liquid. 

The present paper is organized as follows. In § 2 we will calculate the effective potential 
energy felt by the electron using the longitudinal dielectric function of the relativistically de- 
generate electron liquid derived by Jancovici (1962). We will thereby calculate the screening 
potential which will be used for the calculation of the screening corrections to the electron 
capture rates. In § 3 we will summarize the numerical results in the form of analytic fitting 
formulae which will facilitate the application of the present results. We will give concluding 
remarks in § 4. 

2. SCREENING OF THE COULOMB POTENTIAL BY THE 
RELATIVISTICALLY DEGENERATE ELECTRON LIQUID 

The Coulomb corrections to the beta decay rates have been discussed in many standard 
textbooks on beta decay (e.g., Morita 1973). Some authors (Takahashi, El Eid, & Hillebrandt 
1978; Fuller, Fowler, & Newman 1980) have discussed possible importance of the screening 
effects on the electron capture rates at extremely high densities. Therefore, it is extremely 
important to calculate accurately the screening corrections to the electron capture rates by 
the relativistically degenerate electron liquid. 
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We consider the density-temperature regime in which the electrons are strongly degen- 
erate. This condition is expressed as (Itoh et al. 1983) 



5.930 x 10 a 



{l + 1.018(Z/A) 2 / 3 P 6 /3 } 1/2 -l 



[K], 



(2.1) 



where T F is the electron Fermi temperature, p 6 is the mass density in units of 10 6 g cm -3 , and 
Z and A are the atomic number and mass number of the nucleus considered, respectively. 
We assume complete pressure ionization. Therefore, no bound electrons will exist around 
the nucleus. This assumption is valid when the following condition is satisfied (Itoh et al. 
1996): 



-p > 0.378 Z 3 [gem" 3 ], 



(2.2) 



where p is the mass density. 



The potential energy of the electron due to the pure Coulomb potential —Ze 2 /r is 
changed to V(r) because of the screening by the relativistically degenerate electron liquid. 
In this paper we use the static longitudinal dielectric function due to the relativistically 
degenerate electron liquid calculated by Jancovici (1962) based on the relativistic random- 
phase approximation. Jancovici's static longitudinal dielectric function is written in the form 
(Itoh et al. 1984) 
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where n e is the electron number density. The electron Fermi wavenumber is expressed as 



ki 



2.613 x 10 10 T p 6 [cm" 1 



1/3 



(2.8) 



-4- 



The ion-sphere radius Oj is written as (Itoh & Kohyama 1993) 

\™\n { = 1, (2.9) 

1 /3 

cii = 0.7346 x 1(T 10 (Jj) [cm] , (2.10) 
where rii is the ion number density. We also have the relationship 

kpoi = (^) 1/3 ^ 1/3 - (2.H) 

The potential energy of the electron V(r) which takes into account the screening by the 
relativistially degenerate electron liquid is written as 

z ' d 3 k 



2tt 2 J k 2 e(k,0) 

2 r°° , 

- / / \ d Q- (2-12) 

71 Jo " 



Ze 2 (2k F ) 2 r°° sin[(2A; F r)g] 
2k F r n J qe(q,0) 

Therefore the increment of the potential energy due to the screening by the relativistically 
degenerate electron liquid is 

Ze 2 * 



V s (r) = V(r) 



r 



= Ze 2 (2k F )— — <1- - / LV . , HI dq } 

1 F, 2k F r\ 7t7 qe(q,0) H ) 

Equation (2.12) is based on the linear response theory (e.g., Pines & Nozieres 1966). 
Dharma-wardama & Perrot (1982, 1990) have carried out the density-functional study of 
hydrogen plasmas as well as the density-functional study of C, Si, and Ge metallic liquids and 
have found that the detailed results of the density-functional calculations of these systems 
are close to the results obtained by the linear response theory. They have also confirmed that 
the density-functional theory as well as the linear response theory reproduce satisfactorily 
the experimental results on the Ge metallic liquid, thereby proving the excellence of these 
theories for this system. This fact gives great support to our use of the linear response theory 
in the present problem. 

Compared with the case of the pure Coulomb potential, the effective electron en- 
ergy at the position of the nucleus is increased by 1^(0). Hence the usual factor in the 
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electron-capture rates pEF(Z, E) (Morita 1973; Fuller, Fowler, & Newman 1980; Langanke 
& Martinez-Pinedo 2000) should be replaced by 

{ [E - V s (0)] 2 ~ m 2 e } 1/2 [E - V s (0)\ F[Z,E- V s (0)} . (2.14) 

At the same time the shift in the threshold energy for the electron capture should be taken 
into account (Couch & Loumos 1974; Gutierrez et al. 1996; Bravo & Garcia-Senz 1999). 

In Figure 1 we show the function 

t Jo <?£(?. 0) 

as a function of R = 2k pr for various values of r s . In Figure 2 we show the function 

J , 1 U.l f°° MWrr)q] 1 

2k F r\ w J a qe(q,0) v j ' 

as a function of R = 2kpr for various values of r s . In Figure 3 we enlarge Figure 2 for small 
values of R = 2k F r. In Figure 4 we show J(R = 0) as a function of r s . Equation (2.13) can 
be rewritten as 

(Z \ 1/3 

V s (r) = 7.525 x 10~ 3 Z I - P A J [MeV] . (2.17) 

The nuclear radius r nuc can be expressed as (Morita 1973) 

r nuc = l^AV^fm]. (2.18) 

From equations (2.8) and (2.18) we have 

2k F r nuc = 6.3 xlO^Z^p 1 / 3 . (2.19) 

From Figure 3 we see that the function J(R) is almost constant for small values of R which 
corresponds to 2kpr nuc given by equation (2.19), when the mass density is not extremely high. 
This fact justifies our use of Vs(0) for the screening potential at the nuclear radius. When 
the mass density becomes extremely high, then one should use the value of the screening 
potential V s (r) by adopting the value of the function J(R = 2k F r nuc ). 

The electron Fermi energy is given by 

1/2 



Ep = 0.5110 



'l + 1.018(Z/A) 2 / 3 P 6 /3 } - 1 



[MeV] . (2.20) 
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When the electrons are extremely relativistic, we have 

1.460 x 10~ 2 Z J. (2.21) 



Therefore, we find that the screening potential at the origin V s (0) is typically a few percent 
of the electron Fermi energy. Thus we conclude that the screening corrections to the electron 
capture rates by the relativistically degenerate electron liquid are not as great as anticipated 
by Takahashi, El Eid, & Hillebrandt (1978) and also by Fuller, Fowler, & Newman (1980). 

In passing, it is interesting to compare our present detailed calculation with the Fermi- 
Thomas model which is generally more crude than the random-phase approximation that 
has been adopted in the present paper (Pines and Nozieres 1966). In the Fermi-Thomas 
model the screening potential at the origin is expressed as (Pines and Nozieres 1966) 

[V S (0)] FT = Ze 2 k FT , (2.22) 

where kpr is the Fermi-Thomas screening wave number. When the electrons are extremely 
relativistic, the Fermi-Thomas screening wave number is expressed as (Flowers & Itoh 1976) 

kpr 1 74 1 ^ 1/2 



2k F ' 2 



/4 1 

= 0.0482. (2.23) 

\n 137.036; V ; 



By Comparing with Figures 2, 3, 4, we find that the prediction of the Fermi-Thomas model 
is quite good (the accuracy is about 7%) for the present problem when the electrons are 
extremely relativistic. 



3. ANALYTIC FITTING FORMULAE 

In this section we will present accurate analytic fitting formulae for J(R) in order to 

facilitate the application of the numerical results obtained in the present paper. We have 

carried out the numerical calculations of J(R) for 0.00001 < r s < 0.1, < R < 50.0. We 
express the analytic fitting formula by 

10 

J{r s ,R) = J^Oij-sV, (3.1) 

i,j=0 

s = l - (log 10 r s + 3) , (3.2) 

U = 2L {R ~ 25 - 0) • (3 ' 3) 
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The coefficients are presented in Table 1. The accuracy of the fitting is generally better 
than 0.1%. 

The value of J(r s , R) for R = is of course obtained by the fitting formulae (3.1), (3.2), 
(3.3). For the sake of simplicity we will give a separate fitting formula 

10 

J(r s , R = 0) = (3.4) 

i=0 

s = \ ( lo gio r s + 3) • (3.5) 

The coefficients b{ are presented in Table 2. The accuracy of the fitting is generally better 
than 0.1%. 



4. CONCLUDING REMARKS 

We have studied the screening corrections to the electron capture rates by the relativis- 
tically degenerate electron liquid. In particular, we have calculated the screening potential 
caused by the relativistically degenerate electron by using Jancovici's (1962) static longi- 
tudinal dielectric function. We have found that the screening potential is typically a few 
percent of the electron Fermi energy. Hence we conclude that the screening corrections to 
the electron capture rates at high densities are not as great as anticipated by Takahashi, El 
Eid, & Hillebrandt (1978) and also by Fuller, Fowler, & Newman (1980). 

We have presented accurate analytic fitting formulae which will be useful when one 
wishes to apply the present results to the calculations of the screening corrections to the 
electron capture rates at high densities. 

We wish to thank K. Nomoto for drawing our attention to the present important prob- 
lem. We are grateful to Y. Takada for clarifying the validity of the linear response theory 
in the present problem. We also thank Y. Oyanagi for allowing us to use the least-squares 
fitting program SALS. This work is financially supported in part by Grants-in-Aid of the 
Japanese Ministry of Education, Culture, Sports, Science, and Technology under contracts 
13640245, 13740129. 
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Fig. I.— I as a function of R = 2k F r for r s =0.1, 0.01, 0.001, 0.0001, 0.00001. 
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Fig. 2.— J as a function of R = 2k F r for r s =0.1, 0.01, 0.001, 0.0001, 0.00001. 
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Fig. 3. — Same as Fig. 2, but for small values of R = 



